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Abstract. We used the KOSMA 3m telescope to map the core 7' x 5' of the Galactic massive star forming region 
W3Main in the two fine structure lines of atomic carbon and four mid- J transitions of CO and ^^CO. The maps 
are centered on the luminous infrared source IRS 5 for which we obtained ISO/LWS data comprising four high- J 
CO transitions, [Cll], and [Ol] at 63 and 145 fim. In combination with a KAO map of integrated line intensities 
of [C ll] (Howe et al. 1991), this data set allows to study the physical structure of the molecular cloud interface 
regions where the occurence of carbon is believed to change from C^ to C", and to CO. The molecular gas in 
W3 Main is warmed by the far ultraviolet (FUV) field created by more than a dozen OB stars. Detailed modelling 
shows that most of the observed line intensity ratios and absolute intensities are consistent with a clumpy photon 
dominated region (PDR) of a few hundred unresolved clumps per 0.84 pc beam, filling between 3 and 9% of the 
volume, with a typical clump radius of 0.025 pc (2.2"), and typical mass of 0.44 Mq. The high-excitation lines of 
CO stem from a 100 — 200 K layer, as also the [C l] lines. The bulk of the gas mass is however at lower temperatures. 

Key words. PDR - interstellar medium: atoms - interstellar: matter - stars: formation - individual objects: W3 



1. Introduction 

Photon Dominated Regions (PDRs) are surfaces of molec- 
ular clouds where heating and chemistry are dominated 
by FUV (6eV < hv < 13.4 eV) radiation. The photons 
can scatter deep into the molecular clouds due to their 
clumpy structure. A large fr action of the volume and mas s 
is thus photon dominated i Tielens &: Hollenbachlll985|) . 
While FIR continuum emission by dust is by far the ma- 
jor coolant, the gas cools mainly via the atomic fine struc- 
ture lin es of O I, C ii, C i, and the rotational CO lines (e.g. 
iKaufm an et al.lll999j) . Carbon bearing species are thus of 
great importance. The occurrence of carbon with depth 
into the surfaces changes from C"*" to C° and -finally- 
to CO. [Cii], [Ci] and mid- J CO Hues trace the outer 
layers of photon dominated molecular clumps. However, 
[C ll] emission also stems from H li regions, H i clouds (i.e. 
the cold neutral me dium (CNM) (Kulkarni & Heilcs 1987; 
IWolfire et al.ll995|) '). and diffuse H ii r egions (i.e. the warm 
ionized medium (WIM) l)Heileslll99^ '). In the context of 
PDRs, depending on the induced chemistry, the critical 
densities, and energy levels of their transitions, [O i] , [C ii] , 
[C i] , and CO lines serve as probes to trace different tem- 
perature and density regimes. In particular the two fine 



structure lines of neutral atomic carbon at 492 and 809 
GHz provide information on the intermediate regions be- 
tween atomic and molecular gas while tracing two different 
physical regimes. 

Few studies have so far tried to map the major cool- 
ing lines of C ii/C i/CO in individual regions to derive 
a coherent picture of the chemical, physical, and -using 
velocity resolved line profiles- dynamical structure of the 
emitting regions. The studies of the massive star form- 
ing regions SI 06 by ISchne ider ct al.' ("2003', '2002') and 
Orion Molecular Cloud 1 by Mookcriea ct al. ( , 2003^) and 
lYamamoto et alJ l|200lh , as well as the study of the reflec- 
tion nebula NGC 7023 by ICeri n fc PhiUipj l)l998l) show 
that a coherent picture is not easily found using standard, 
steady-state PDR models. Especially the large spatial ex- 
tension of neutral carbon and the high abundance deep in- 
side molecular clouds appears to be in contradiction with 
standard PDR m odels. This is also a ddressed by the [C l1 
vs. CO studies of IPhime et a l.' ('1994') in S 140.'O kaet al 
1I2OOII) in DR15.lKameg-ai et al. (2003') in p Ophiuchi, and 
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lOka et al.1 l)2004D in NGC 1333. In a study of t he tra nslu- 
cent cloud MOLD 123.5 + 24.9, iBensch et aP l)2003 l) find 
lower than expected [C i] int ensities and JC l] t o CO ratios. 
iKameeai et"al] ()2003|) and lOka et all l)2004|) argue that 
[C i] traces young evolutionary stages of cloud surfaces 
which is consistent with time-dependent model calcula- 
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Table 1. List of observational parameters. The columns give line frequency, main beam efRcieny i?eff , half power beam 
width, obsering mode, number of mapped positions, velocity channel width Av (the resolution is a factor of ^ 1.4 
larger) , mean atmospheric zenith opacity (tq*™), and observing period. The forward efficiency F^s from skydips is 
constant at 0.9. 



Species 


Transition 


Frequency 

[GHz] 


Soft 


HPBW^ 


Mode 


Pos. 


Av 
[m s-i] 


{rS'n 


Observing period 


i^CO 


J = 3-2 


345.7960 


0.65 


76 


OTF" 


1091 


29 


0.54 


11/2000 


"CO 


J = 3-2 


330.5880 


0.65 


76 


OTF 


1091 


29 


0.63 


1/2001, 8/2002 


i^CO 


J = 4-3 


461.0408 


0.57 


57 


DBS'' 


192 


675 


1.09 


10/2001 


i^CO 


J = 7-6 


806.6518 


0.25 


42 


DBS 


192 


389 


1.15 


12/2001 


"CO 


J = 8-7 


881.2728 


0.38'* 


40'* 


DBS 


60 


356 


0.89 


1/2003 


12 [C I] 


^Pi-^Po 


492.1607 


0.50 


55 


DBS 


192 


634 


0.82 


12/2001 


12 [C I] 


^P2-^Pi 


809.3420 


0.25 


42 


DBS 


192 


388 


1.10 


12/2001 



OTF - On-The-Fly observing mode using an emission-free off position at (42.5', —25.7') 
DBS - Dual-Beam-Switch mode using the secondary mirror (6' throw in azimuth) 
The main beam efficiency and the HPBW were determined from cross scans of Jupiter. 

Data at 880 GHz were taken after an adjustment of the primary mirror following photogrammetric measurements by 
iMiller et al.l i2002f) . The beam efSciency at 880 GHz was extrapolated from new efficiencies taken at 810, 492, and 345 GHz 
using the Ruze formula, consistent with a surface accuracy of 25 fim rms. The HPBW was scaled using the width at 810 GHz. 



tions (cf. also lStorzer et alJl997|) . lGerin fc PhillipsI l)l998|) 
suggest that a second PDR deeper in the cloud explains 
the observed enhanced [Ci]/[CO] ratios in NGC7023. The 
importance of the interclump medium, providing effective 
shi elding of the CO in the embedded clumps, is discussed 
bv lBenschet all (|2003|) . 

To improve our understanding of PDRs in different 
sources with different geometries and different FUV radi- 
ation fields, we present here a study of W3 Main. As the 
major cooling lines with their bright transitions are partic- 
ular important as tracers in external galaxies, proper and 
detailed understanding of their emission and the chemical 
and physical conditions of the originating gas is necessary 
as a firm basis for their interpretation. 

The chain of giant Hii regions W3-W4-W5 lie at 
the edge of the Cas 0B6 association and for m part of 
the Perseus spiral arm. lHever &: TerebevI (|1998|) combined 
large-scale CO 1-0 FCRAO maps wi th Hi line emission 
and 21cm continuum DRAO data (Nor mandeau et alJ 
Il997|) fcf. lRead 1981) . and IRAS HIRES images at sim- 
milar angular resolutions of about 1'. W3 is associated 
with the W3 giant molecular cloud compri sing more than 
lOfM© at a distance of 2.3 kpc (Georgelin &: Georgeliri 
Il976l) . This distance wi ll be u sed throughout this pa- 
per, though llmai et alJ (2000; recently found indica- 
tions for a distance of only 1.8 kpc. W3Main is a well 
known site of hig h-mass star formation comp r ising sey - 
eral Hii regions JTieftrunk et all 1199,4 Il997t ll998b 'aV 
ISO/SWS an d LWS spectroscop y of the compact Hii 
region W3A l|Peeters et al.l l2002|) reveals strong atomic 
emission lines. W3A is associated with a luminous in- 
frared source IRS 2, belon ging to a cluster of i nfrared 
sources (IRS 1-7) found by Iwvnn-Wimams et alJ 1^72); 
iRichardson et all l)l989(l . The total luminosity of IRS 5 
is estimated to be 2 - 5 10^ Lq (|Werner et all Il980l: 



iThronson et al.lll980t fCampbell et"al]ll995|) . [Ladd et alJ 
l|l993^ observed the associated emission of dust in the 
submillimete r identifying three cor es (SMS 1—3). NIR ob- 
servations bv lMegeath et al.l 1)199611 reveal that high-mass 
star formation is associated with the formation of a dense 
cluster of several dozens of low-mass stars. 

Column de nsities of molecular gas were mapped by 
iTieftrunk e~l. (1995) in mil limetric C^^O and C^^S hues 
and bv iRoberts et alJ ^ll997^ in ^^CO and Ci^O 1 - 0. 
The latter observation were conducted with BIMA and 
reveal six clumps with diameters between 0.1 and 0.3 pc 
(masses between 120 and 480 M©) in the vicinity of IRS 5 
and IRS 4. The denser and warmer molecular g as in the 
vicinity of IRS 5— IRS 4 was observed bv Hasegawa et alJ 
(1994) in CO 3 - 2 and 6 - 5. 

As part of a larger stud y of dense, hot g as in 
Galactic star forming regions, iKrtigel et alJ l)l989tl and 
iBoreiko fc Bet3 lll99ll) observed ^^CO 7 - 6 and ^^CO, 
^^CO 9 — 8 towards W3IRS 5, respectively. The large-scale 
extent of the 492 GHz line of ato mic carbon an d of ^^ CO 
and C^^O 2 - 1 was observed bv Flume et alJ l|l999l) at 
3' resolution. How e et al.l l|199lf) mapped the [C ll] 158 fim 
emission of W3 Main and developed models to reconcile 
the distribution of [C ii] emission and the FUV field. 

Here, we attempt to model the mapped emission of 
mid- J and high- J CO lines, of both [C i] lines, and of [C ii] , 
at resolutions of ~ 1'. Section 2 describes the KOSMA ob- 
servations and the ISO/LWS data set. Section 3 presents 
the maps of integrated intensities, spectra at four selected 
positions, and a first derivation of physical properties. In 
Sect. 4 we discuss the physical structure of the photon 
dominated regions. We present a model of the FUV field, 
needed for a detailed comparison of line integrated intensi- 
ties with PDR models. We summarize our results in Sect. 
5. 
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Fig. 1. (Left) Integrated ^^CO and (Right) ^^CO 3—2 large scale maps of W3, including WSMain in the north and 
W3(0H) in the south. The integration interval is [—55,-25 kms~^] and thus also covers the emission of W3(0H). 
Offsets are given in arcsec. Contour levels are described in the upper left corner of each map: first contour= 3 cr/step 
size/peak intensity. The position (0, 0) of IRS 5 is marked by a star, IRS 4 is marked by a square. Triangles mark the 
positions of IRS 1 (to the cast), IRS 3, and W3(0II). The inner black box indicates the size of the higher frequency 
maps shown below. The 0/0 position is at IRS 5: a(1950) = 02^21^53!2 (5(1950) = 61°52'21". 
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Fig. 2. (From upper left to lower right) Integrated intensity maps of CO 4—3, CO 7—6, [C i] 1-0 and [C i] 2—1 of the 
W3 Main region (indicated by the rectangle in Figure^). 
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Fig. 3. Map of "CO 8-7 spectra at a resolution of 75" 
(0.84 pc) and on the Tmb scale after a pointing correction. 
The temperature range is —2 to 10 K. The visi range is 
-70 to -lOkms-i. 



2. Observations 

The KOSM A 3m submillimete r telesc ope on Go rnergrat, 
Switz erland l|Winnewisser et al. 1986: Kramer et al. 1998. 
I2OO0I) was used to observe the transitions of carbon and 
the mid- J transitions of CO. 



2.1. [Ci] and mid- J CO 

We have mapped W3 (Fig.^ EJ in the two atomic car- 
bon fine structure transitions ([Ci]) at 492 GHz (609 /zm, 
3Pi-3Po; hereafter 1-0) and 809 GHz (370 ^im, ^Pa-^Pi; 
hereafter 2—1) and the mid- J rotational transitions J — 
4-3 and J = 7-6 of CO, and J = 8-7 of "CO using the 
newly installed SubMillimeter Array Receiver for Two fre- 
quencies (SMART) (Gra f et al. 2002^ on KOSMA. 

SMART is a dual-frequency eight-pixel SIS-heterodyne 
receiver that simultaneously observes in the 650 and 
350 /im atmospheric windows (jGraf et al.ll2002.'l . It con- 
sists of two 4 pixel subarrays, one operating in the fre- 
quency band of 455 — 495 GHz and the other in the band 
of 795 - 882 GHz. The IF signals are analyzed with two 4- 
channel array-acousto-opti cal spectrometers with a spec- 
tral resolution of 1.5 MHz iJHorn et al.l 'l999'l. The IF fre- 
quency of 1.5 GHz is ideally suited for simultaneous detec- 
tions of the [Ci] 2-1 line at 810 GHz and the CO 7-6 line 
at 806 GHz in opposite sidebands. The two subarrays are 
pointed at the same positions on the sky. Thus, four spa- 
tial positions are observed simultaneously. The [C i] 1-0 & 
2—1 and the CO 7—6 observations of W3 presented here 
were done simultaneously. The relative pointing is thus 
constant. Due to the simultaneous observations at both 
frequency bands, relative calibration errors are strongly 
reduced and line r atios are measured precisely (cf. e.g. 
IStutzki et alJll997(l . The typical receiver noise tempera- 
tures achieved at the center of the bandpass is around 
150 K at 492 GHz and between 500 and 600 K at 810 GHz. 
The receiver setup consists of a K-mirror type image ro- 
tator, two Martin-Puplett diplexers, and two solid state 
local oscillators, which are niultip lexed using coUimating 
Fourier gratings JGraf et alJ2002|) . The image rotator cor- 
rects for the sky rotation during the observation and it 
allows to set the array to any desired orientation. 



The observations were performed in Dual-Beam- 
Switch mode (DBS), with a chop throw of 6' in azimuth. 
This mode of observation produces flat baselines. Only a 
zeroth order baseline was subtracted from the calibrated 
spectra. The spacing between two adjacent pixels is 114". 
Pointing was frequently checked on Jupiter and the sun. 
The resulting accuracy is better than ±20". 

An area of ^ 7' x 5' was mapped in the two [C i] lines 
and the two CO lines, on a 27.5" x 27.5" grid with a total 
integration time of 160 s per position. Due to the extent of 
the W3 cloud complex weak self-chopping artefacts were 
found at the borders of the mapped region. This was con- 
firmed by separate observations of these positions in total 
power position switch mode. We did not correct for these 
artefacts. 

Atmospheric calibration was done by measuring the 
atmosph eric emission a t the OFF-position to derive the 
opacity l|Hivanial Il998|) . Spectra of the two frequency 
bands near 492 and 810 GHz were calibrated separately. 
Sideband imbalanc es were corrected using standard at- 
mospheric models (|Cernicharcilll985|) . The CO 7—6 line 
observed in the image sideband, was also corrected. 

A smaU map of the "CO 8-7 line at 880 GHz covers 
an area of 5.5' x 2.3' centered on IRS 5 and IRS 4 (Fig.EJ. 
Integration times are on average 480 sec per position. No 
self-chopping effects are discernable in the resulting spec- 
tra. A pointing offset was corrected off-line by shifting 
the "CO data by (0",-27.5"). The main beam efficiency 
BcS and HPBWs were deter mined using c ontinu um cross 
scans on Jupiter (Tabled). ICriffin et all l|l986l Table 1) 
observed Jupiter and Mars at 10 frequencies in the 0.35 
to 3.3 mm wavelength range using Mars as primary stan- 
dard. We used their Jovian brightness temperatures, lin- 
earily interpolated, to derive the temperatures at the fre- 
quencies of the KOSMA observations. The KOSMA half 
power beamwidths (HPBWs) were derived from the ob- 
served FWHMs. At 810 GHz - 25% of the emission is 
detected by the main beam. Continuum scans of the edge 
and disc of the sun reveal an extended Gaussian error 
beam component of -- 200" FWHM, detecting about 40% 
of the emission. We attribute this error beam to resid- 
ual misalignment of the 18 panels of the primary mirror. 
Due to the large extent of the errorbeam relative to the 
mapped region, we cannot convolve out the errorbeam. 
Instead, we simply scaled the antenna temperature data 
with the ratio of forward over main beam efficiency to go 
from the T^ scale to the T^^h scale. 

From observations of reference sources, we estimate the 
relative calibration error to be ^ 15%. 



2.2. Low- J CO & ^^CO 



We have used a dual-channel SIS receiver IjGraf et alJ 
1 19981') on KOSMA to observe large scale maps in CO 
J = 3-2 and "CO 3-2 emission of W3 l|Miilleil 12001]) . 
The obser vations were done using th e On-The-Fly ( OTF) 
technique (JBeuther et al. 2000: Kramer et al.l 
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fully sampled grid of 30' x 30'. For the low- frequency obser- 
vations we used the high resolution spectrometer (HRSl) 
till Winter 2001 and the variable resolution spectrometer 
(VRSl) lateron. Sinusoidal baselines were subtracted. See 
Tabled for a summary of observational parameters. 



2.3. ISO/LWS spectroscopic data 



ISO/LWS l)Clegget alJll996|) spectra are available from 
the data archive for the position of IRS 5 in W3 Main 
(TDT 47301305). The grating spectrometer covers the 
wavelength range 47— 197/xm. The standard OLP-10 data 
product was used along with ISAP v. 2.1 to deglitch, de- 
fringe, and average the data. Linear baselines were sub- 
tracted from all spectra. Gaussian fits, with the line width 
held fixed at the instrumental width, were done to derive 
integrated intensities. 



3. Results 

3.1. Why to use different tracers? 

The observed transitions trace varying physical condi- 
tions of the emitting regions. Transitions are characterized 
by their critical density for collisions with H2 or other 
partners and by their upper level energy (see Table 2 in 
iKaufman et alJll9 99). For the CO J ^ J — 1 transition, 
the critical density is 4 10^ J^ cm~^ and the upper level 
energy corresponds to 2.8 J{J+1) Kelvin. The mid- J CO 
transitions observed with KOSMA are thus sensitive to 
gas at about 10^ to 10^ cm~^ at temperatures between 30 
and 200 K while the high- J hues traced by ISO/LWS are 
sensitive to typically lO^cm"^ and 650 K. The regime of 
high densities > 10^ cm~^ and high temperatures > 200 K 
is also traced by the two atomic [O i] lines at 145 /xm and 
63 /zm. On the other hand, the two atomic [C i] lines and 
the [Cii] line trace energies of less than 100 K and densi- 
ties of less than 3 10^ cm~'^. 



3.2. Maps of integrated intensities 

We present here KOSMA maps of W3 Main at spatial res- 
olutions of about one arcminute. The ^^CO and ^'^CO 3 — 2 
maps include W3Main and W3 (OH), and cover an area 
of ~ 20' X 14'. A 7' X 5' region centered in W3Main-IRS5 
was mapped in [Ci] ^Pi-SPo, ^Pg-^Pi, ^^CO 4 - 3, and 
7-6. The 5.5' x 2.3' map of ^^CO 8-7 covers the IRS 5 
and IRS 4 posi tions. We have i ncluded in our analysis the 
[C 11] map of H owe et al.l l)l99l|) and ISO/LWS data at the 
position of IRS 5. 

In contrast to the velocity integrated CO maps which 
show only one maximum near IRS 5, the corresponding 
spectra show many details like self-absorption features and 
overlapping velocity components. These will be discussed 
in Sect.rOl 



3.2.1. Maps of 12C0 and ^'KO 3-2 

The ^^CO and ^'^CO 3—2 maps exhibit («) a single central 
peak close to the position of IRS 5, (ii) an extension t o 
the north-east toward W3 North fcf. IPlume et alJll999l) . 
which is the northernmost member of the string of active 
star-forming regions, and (iii) a ridge of emission extend- 
ing towards the south in the direction of W3(0H). In 
addition, (iv) both maps show another ridge of emission 
extending eastward from IRS 5 and then bending north- 
east. The two north-eastern ridges are more clearly seen 
in the ^^CO isotopomer than in ^^CO. 



3.2.2. Maps of CO 4 - 3 and 7 - 6 

The higher-frequency maps cover the central part of 
W3Main encompassing an area of ~ 7' x 5'. 

The CO 4 — 3 and 7 — 6 maps again show a strong peak 
near IRS 5, slightly shifted towards IRS 4, which drops 
rather steadily in all directions. A ridge of emission ex- 
tending from IRS 4 in south-western direction is discern- 
able in the CO 7 — 6 data set. 



3.2.3. Maps of [Ci] 1 - and 2 - 1 

Comparison of the C i maps with the CO maps in Figures^ 
and [3 shows that although there is a large-scale similarity 
between the two emission patterns, the mid-J CO emis- 
sion is more centrally concentrated than the [C i] emission. 
This indicates the origin of the CO emission in regions 
with high density and temperature, while the [C l] emis- 
sion originates in regions of comparatively lower density 
and temperature, and hence traces out the widespread dif- 
fuse UV illuminated neutral gas. 



3.3. Spectra at selected positions 

In FigureQI we present spectra at four positions, IRS 4, 
IRS 5, and two positions 0.92 pc (85") to the north and 
to the south of IRS 5 where the line integrated intensities 
drop to about a third of the peak intensities. In general, 
^^CO line shapes show a complicated structure, espcially 
at IRS 5 and 4, while ^'^CO and [C i] line shapes show little 
deviations from Gaussian profiles. 

At the position of IRS 5, the ^^CO 3-2 and 4-3 spectra 
clearly show an absorption dip, while the ^^CO 7-6 spec- 
trum is fiat-topped. The ^^CO and [Ci] spectra all peak 
at about the velocity of the dip, near —38 kms~^. The 
foreground absorption in the center region of W3 Main 
is alread v well known from the CO observations of 
IHasegawa et al.l l|l994l) . It is also seen in the water line 
observed with the Odin satellite l|Wilson et al.ll20fl3l) . In 
order to correct for self-absorption in the ^^CO line pro- 
files, we have fitted Gaussian profiles to the line wings 
(Fig. 21 to derive reconstructed integrated intensities and 
their ratios (Table[21). In the following, we will only use 
the corrected line integrated ^^CO intensities. 
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Fig. 4. Spectra at four representative positions, including IRS 5 and IRS 4. All spectra are on a common spatial 
resolution of 75" (0.84 pc) and are plotted on the Tmb scale at 0.1 kms^^ channel spacing. ^"^CO and [Ci] lines peak 
at the velocity of a dip in the ^^CO spectra, indicating that the dip is due to a colder foreground component at about 
—38 kms~^ (denoted by a vertical line). To estimate the line integrated ^^CO intensity of the background source, we 
fitted a Gaussian to the line wings of the ^^CO spectra. FigureEl gives an example. 



The i^CO FWHM line widths at IRS 5 are - 9 kms-\ 
larger than those of the "'^'^CO lines which are in turn larger 
than those of the [Ci] lines which exhibit ^ 6 kms^^. 

3.4. Temperatures and carbon budget at four 
representative positions 

Observed line integrated intensities at the same angular 
resolution of 75" and some derived line ratios are compiled 
in Tables|21and|31 Physical parameters at the four selected 
positions are listed in TableQ] 



3.4.1. CO line ratios, excitation temperatures, and 
beamfilling factors 

The i^CO 7~6/i2co 4-3 line ratio of integrated inten- 
sities (TableE)), corrected for self-absorbed hue profiles, 
varies between 1.0 and 1.3. This is consistent with opti- 
cally thick emission. Assuming equal excitation temper- 
atures (LTE), a ratio of 1±20% indicates temperatures 
of more than 40 K. This is consistent within the errors 
with the excitation temperatures derived from ""^^CO 3- 
2 peak line temperatures in the center region and to the 
North, again assuming optically thick emission (Table|2|). 
Excitation temperatures from CO 7-6 are ~ 10 K higher 
due to the Rayleigh- Jeans correction at equal peak line 
temperatures. 

Only the ISO/LWS observations of high- J CO transi- 
tions at IRS 5 together with an escape probability analy- 
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Fig. 5. The ^^CO 3-2 spectrum at IRS 5 is shown to- 
gether with the Gaussian line profile fitted to the line 
wings. The ratio of fitted to observed line integrated inten- 
sity is 1.5. The lower box shows the residual of the profile 
fit. 



sis ('ChaDter l4.4.3(l show that an additional gas component 
exists exhibiting temperatures of more than 140 K. 
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These arguments, assume so far tfiat radiation is emit- 
ted from a homogeneous region. However, temperature 
and density gradients are known to exist in UV iUumi- 
nated clouds. We wiU explore their significance when com- 
paring the observations with PDR models in Sect. 14. 51 

3.4.2. [Ci] line ratios 

Ratios of the two [C i] lines vary systematically between 
1.3 in the south and 2.9 in the north (Table|2J). These 
ratios indicate optically thin emission since for optically 
thick [C i] lines we would expect line ratios of about 1 
or less. The low peak line temperatures of less than 11 K 
do not reflect the kinetic temperatures estimated above, 
instead they also indicate low optical depths (t < 0.15 for 
Tex = 100 K). 

In the optically thin limit, the excitation temperature 
is a sensitive function of the line ratio R in Kkms^^: 
Tex = 38.8K/ln[2.11/i?]. In Column (8) of Tabled we 
give the lower limit on T^x using this formula given the 
observational error on the line ratio. The error is such 
that only lower limits on excitation temperatures can be 
derived. Ratios above 2.11 as found at one position north 
of IRS 5 lead to unphysical negati ve excitation ternpera- 
tures. Detailed PDR modelling bv lKaufman et alJ ijlQQIJ . 
Fig. 8) confirms that for densities of upto n = 10^ cm~'^ 
and FUV fields of upto Go = 10^-^ the CI line ratio stays 
below 2.1. 

To exclude systematic observational errors, e.g. due to 
the KOSMA error beam at short wavelengths or due to 
the atmospheric calibration, we studied the distribution 
of all observed CI line ratios. In total 167 positions were 
observed, out of which 14 positions (8.4%) show a CI hne 
ratio greater than 2. This is consistent with a true ratio 
of 1. 5. a 20% Ig error, and a Gaussian error distribution 
fe.g. lTavloiil982ll . Using the above simplifying formula for 
optically thin emission and equal Tox, a ratio of 1.5 results 
in a temperature of 114 K which is characteristic for the 
average temperature of the CI emitting gas. 

In conclusion, carbon excitation temperatures derived 
from the [C i] line ratios are typically 100 K, and are higher 
than those derived from ^^CO. However, the sensitivity of 
the above ansatz on the observational errors for high ob- 
served ratios does not allow more quantitative conclusions. 

Only few observations of the upper [C i] line are re- 
po rted so far (T ablcEll. The 809 GHz line was first detected 
by Jaffe et al.l (|li)85) and later observed bv lCenzel et alJ 
l)l988(l in M 17 and W 51. The COBE satellite observed 
the Milky Way at 8° re solution and found low line ra- 
tios (jFixsen et al.l Il999l) of 0.31 in the inner Galactic 
disk, corresponding to a kinetic temper atures of > 20 K. 
Line r atios observed in Orion A ('OMCl) lYamamoto et al.l 
l|200l|) using the Mt.Fuji telescope indicate excitation tem- 
peratures covering the range of 52 to 62 K and more. 
IZmuidzinas et al.l l)l988() detected the upper [Ci] line in 
several Galactic star forming regions with the KAO, and 
derived excitation temperatures between 25 K and 77 K 



from the [Ci] 2-1/1-0 line ratios. They derived a carbon 
excitation temperature of 25 K from observations of W3. 
This is in contradition to our result which is however based 
on simultaneous observations with the same instrument. 
Moreove data were smoothed to the same spatial resolu- 
tion thus reducing possible calibration errors. 

3.5. The case for hot gas at IRS 5 

iHelmich et al.l (jl99 4^ found rotational temperatures of 
SO2, CH3OH, and H2CO ranging between 64 and 187K 
when using the J CMT with high spatial resolutions. 
Ivan der Tak et al.l l|2000|) find rotational temperatures of 
CH3OH and H2CO of 73 and 78 K. 

High- J lines of CO and ^^CO w = - 1 (4.7 ^m) ab- 
sorption measurements towards IRS 5 with a subarcsecond 
sized beam giv en by the size of th e IR continuum back- 
ground source iJMitchell et al.lll99d at the CFHT) reveal 
a hot gas component at 577 K comprising ~ 50% of the 
total gas column density per beam at densities exceeding 
10^ cm^'^. A colder component at r^ 66 K comprises the 
other half. 

Sohd CO is also detected towards IRS 5 
llSandford et alJ Il988t iLacv et all Il984t iGeballel Il986ft . 
but con tributes only 0.2% to the total CO column 
density l)Mitchell et'aLlll99Cll). ISO-SWS observ ations of 
the ^^C02 ice band profile l)Boogert et al.ll200(l) towards 
IRS 5 show indications of heated ice, consistent with the 
above results. 

Ammonia observations at resolutions of 40" were used 
bv lTieftrunk et all l)l998bl) to derive kinetic temperatures 
at the positions of IRS 5 and 4. While the kinetic tem- 
perature at IRS 4 is estimated to be 48 K in excellent 
agreement with the excitation temperature we have de- 
rived from ^^CO 3-2, the kinetic temperature at IRS 5 
was estimated to be 150 K. 

In conclusion, observations indicate a strong tempera- 
ture gradient peaking at IRS 5 at upto 600 K. The hot gas 
component has a low beam filling factor while the colder 
gas component is more widespread. Variations of optical 
depths and chemistry are also at play. 

3.5.1. Relative carbon abundances 

Table01 lists the relative abundances of the major carbon 
bearing species in the gas phase: Cii, Ci, and CO. These 
are derived from column densities assuming optical thin 
"CO 3-2, [Ci], and [Cii] emission and LTE. For the Cl 
column density we assumed an excitation temperature of 
100 K at all positions. Note that the column density is 
not temperature sensitive: an increase of Tex by a factor 
of 2 only results in a 10% increase of the column density. 
C II column densities are independent of the excitation 
temperature in the thermalized l imit and for te mperatures 
much higher than 91 K (see e.g. ISchneider et al. 2003). 

The total H2 column density is derived from the ^^CO 
3-2 integrated intensity via the CO column density as- 
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Table 2. Observed integrated intensities in Kkms^^ at the four positions in WSMain shown in Figure^ The ^^CO 
line intensities were corrected for self-absorption. For this reason, some of the integrated intensities are identical. 
The error on the integrated intensities observed with KOSMA is estimated to be ^ 15%. [Cii] intensities are from 
iHoweet al1l)l99ll) . 



Aa/A5 



'CO 



^CO 



[Ci] 



[Ciil 



[","] 


3-2 


4-3 


7-6 


3-2 


8-7 


"Pi-"Po 


^P2-^Pl 


P3/2~ Pl/2 


0,0 (IRS 5) 


408.4 


408.4 


432.4 


148.7 


52.9 


38.3 


73.5 


311.3 


-55.0/27.5 (IRS 4) 


408.4 


408.4 


432.4 


158.6 


37.4 


42.5 


78.2 


267.2 


0,82.5 (North) 


231.2 


231.2 


231.2 


63.7 


- 


15.5 


44.3 


148.8 


0, -82.5 (South) 


129.9 


99.4 


132.7 


60.2 


- 


21.3 


27.4 


311.3 



Table 3. Selected ratios of observed integrated intensities in Kkms ^ at the four positions in W3Main which are 
shown in Figure^l The error on the integrated intensity ratios ob served with KO S MA i s estimated t o be ^ 20%. 

'' 'Wilson et alJ 
COBE), f 



For co mparision. wc also list ratios observed in other sources by ° Schneider et al. f200?, KOSMA 
ll2n0l[ HHT) llzha ng ct al. ( 2001, AST/RO) ^' P etitpas & Wilson ( 1998, JCMT), " Fixscn ctahl (|1_ 
IStiitzki et all lll997l HHT). g llsrael fc BaasI ll2002l JCMT, Fig. 2). ^'lYamamoto et all ll200lL Mt.Fuji) 
ll200(i HHT/IRAM 30m), ^ IWeiss et al.l ll2003i IR AM PdBI), "" iTieftrunk et all ll200ll AST/RD). 
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[Ciil 



[","] 


"CO 3-2 


CO 4-3 


"CO 8-7 


[Ci] 1-0 


CO 4-3 
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CO 4-3 


W3 Main: 
















0,0 


2.75 


1.06 


2.81 


1.92 


0.09 


8.14 


1.12 


-55,28 


2.57 


1.06 


4.24 


1.84 


0.10 


6.29 


0.82 


0,82.5 


3.63 


1.00 


- 


(2.86) 


0.07 


9.60 


1.11 


0,-82.5 


2.16 


1.33 


- 


1.29 


0.21 


14.62 


3.27 


S106" 


3.0 


1.03-1.57 




1.0-1.8 


0.12-0.38 


6.4 


0.8 


Orion A 




1.2-3.5'' 




0.83*^ 








Carina MC 










0.14-0.45 






G333.0-0.4, NGC6334A, 


G351.6-1.3 






0.48-0.75™ 








Inner Galaxy*^ 










0.31 


3.54 


4.29 


Galactic Center*^ 




0.11 






0.22 


1.38 


0.71 
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M82 Center 




0.28-0.36' 




1.0^ 








Sample of 15 galaxies® 










0.1-1.2 






Cloverleaf 








o.si*-- 









suming a ^^CO/^^CO abundance of 65 iJLanger fc Penziaa 
Il99(ll and a H2/CO abundance of 8 10"^ 

In WSMain, the relative abundances of Ci/CO are 
rather constant: the ratio is ^ 0.11 at all four positions. 
Similar ratios of ~ 0.1 are found by IZm uidzina s et alJ 
( 1988) in a s t udy o f seven Galactic star forming regions. 
IPlume et al.l l|2000|) derive similar CI/CO ratios in the 
Orion A cloud, but observed an increase towards the cloud 
edges to values of ^ 0.5. Schneider ct al. (2003) find low 
ratios in S 1 06 var ying between 0.09 and down to 0.03. 
Howe et al.' NOO") find ratios of about 0.4 in M17SW. 



Gerin fc Phillips (1998) find ratios of between 0.1 and 0.7 



in the photodissociation region associated with the reflec- 
tion nebula NGC 7023. The quiescent translucent and thus 
less well shielded cloud MCLD 123.5 +24.9 exhibit s high 
C i/CO ratios of between 0.2 and 1.1 iJBensch et alJl2003i) . 
Back to W3 Main: The C i abundance relative to both 
other two major gas-phase carbon species is constant at 



~ 8%. The Cii abundances is high in the south (32%) 
and lower at the other positions. In other words, the 
relative CO abundance is greater than 60% through- 
out. Similar fractional abundances of Cii:Cl:CO as a t 
W3 Main-South are seen in IC 63 l|Jansen et al.l Il996|) : 
37:7:56, stil l lower abundances of CO are see n in NGC 2024 
(Orion B) iTaffe fc Plumelll995l) : 40:10:50. iBrooks et alJ 
l)2003J) found fractional abundances for Cil:Ci:CO of 
68:15:16 in jyCarina. On the other hand, CO a bundances 
are much higher in S106 ijSchneider et alJl2003|) at all po- 
sitions studied (> 86%), dwarfing the contributions from 
Cii and Ci. In S106, Cii contributes 8% to the carbon 
content at the position of the exciting star S106IR, even 
less at all other positions. In contrast, the Cii fractional 
abundance is higher all over W3 Main, probably due to 
the much larger number of embedded OB stars. 
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Table 4. Physical parameters at four positions in WSMain. The excitation temperature Tex(CO) (col. (1)) is derived 
from ^^CO 3-2 peak temperature. Lower limits of the Ci excitation temperature (col. (8)) are derived from the [Ci] 
2-1/1-0 line ratio assuming a calibration error of 20%. For the calculation of CI column densities, it was however 
assumed that Tex([Ci]) is lOOK constant. Column densities per 75" beam (0.84 pc) and the relative abundances of 
the three major gas phase species CO, Ci, and Cll, are derived assuming LTE (cf. chapter [3. 5. l|l . The FUV field Gq 
smoothed to 75" is estimated from the luminosity of the embedded OB stars (cf. chapter [4. l|l . 



Aa/A<5 


Tex (CO) 

[K] 


iV(CO) 


iV(Ci) 
[10^' cm-^ 


iV(Cn) 


Cl:CO 


Cn:Cl:CO 


Tex (C I- Ratio) 

[K] 


A^(H2) 
[10^1 cm-2] 


Go 
[10*] 


0,0 


53. 


55.2 


5.8 


13.8 


10:90 


18: 8:74 




> 122 


39.9 


40.4 


-55.0,27.5 


53. 


58.7 


6.4 


11.8 


10:90 


15 


8 


76 




> 108 


42.5 


19.9 


0,82.5 


39. 


20.9 


2.4 


6.6 


10:90 


22 


8 


70 




- 


15.1 


10.4 


0,-82.5 


21. 


26.2 


3.2 


13.8 


11:89 


32 


7 


61 




>54 


18.9 


6.7 



3.6. Total H2 column densities and masses 

Total H2 column densities per 75" beam derived from 
""^■^CO vary between 1.5 and 4.3 10^^ cm~^ at the four posi- 
tions (TableQl). This agrees within a factor of 2 with total 
column densitie s derive d from C^^O at 40" resolution by 
iTieftrunk et alJ l|l998h(l . 

The total mass per 75" beam from ^■^CO at IRS 5 is 
481 Mq. Here, we assume a molecular mass per hydrogen 
atom of 1.36. 



4. The physical structure of W3 Main 

4.1. Modelling the spatial structure of the FUV field 

ITieftrunk et alJ l|1997f) presented a multi-frequency VLA 
study of the radio contiuum sources in W3 Main. This re- 
gion exhibits the diffuse regions W3 H, J, and K, the com- 
pact regions W3A, B, and D, the ultracompact regions 
W3 C, E, F, and G, and hypercompact regions with diam- 
eters of less than 1" (~ 0.01 pc) W3Ma-g and Ca which 
probably correspond to the youngest evolutionary stage 
of an H II region. FigureEl shows the po sitions of the 17 
OB s tars exciting the H 11 regions fTablel^T Tieftrunk et alJ 
Il997|) . The most luminous source is IRS 2, identified as an 
05 star. The luminosity o f IRS 5 corresponds to a clus- 
ter of ~ 7 BO- type stars dClausscn ct al.lll993) . We use 
the spectral types of the exciting OB stars derived by 
iTieftr^mk et alJ ( 199711 to calculate the spatial large-scale 
variation of the FUV field heating the molecular cloud. 

The FUV field of W3 Main was modelled in two di- 
mensions by summing the contributions from all sources 
at each pixel of a regular grid. For all sources, we assume 
that the positions inside the cloud are given by their pro- 
jected distances. The FUV field was then derived assuming 
that the stars are ZAMS stars at the e ffectiv e temperature 
corresponding to their spectral type l)Panagiail973) . We 
furthcron assume only geometrical dilution of the flux, i.e. 
we neglect scattering and absorption of UV by dust and 
clumps, thus deriving upper limits only. 

Due to the many bright embedded sources, the result- 
ing FUV field, smoothed to 75" resolution, varies by a 
factor of 100, dropping from Go = 10^'^ in the center to 



lO'^-^ at the edges of the mapped region. The Habing-field 
Go used throughout this paper is in units of 1.6 10^'^ erg 
s-icm-2 (Habine 1968). Sec Drainc & Bertoldi (1996) for 
a discussion of other measures of the average interstellar 
radiation field. 

Th e total FIR luminositv of W3 given in the literature 
e.g. bv lThronson et alJ l)l980|) is 1 10^ Lq. Solely the lumi- 
nosity of the 05 star powering IRS 2 accounts already for 
68% of this (cf. TableEl). The total luminosity of aU OB 
stars is 1.5 10® L©. This indicates that the primary source 
of the FIR continuum is the reemitted stellar radiation 
field of the OB stars and that some of the stellar photons 
escape without impinging on the clouds. Additional heat- 
ing sources like e.g. outflow shocks, stellar winds, or the 
embedded low mass objects are not needed to explain the 
observed total FIR luminosity. 

4.2. A map of [Cii] emission 

FigureEl sho ws the extende d [Cii] ^P3/2— ^Pi/2 158 fim 
emission ftHowe et al.l Il99lh in contours at a resolution 
of 80". The [Cii] emitting region shows a north-south 
extension and extends over all the H 11 regions which 
have been identified in radio continuum emission studies 
IjTieftrunk et al.lll997r . Interestingly, the north-south ori- 
entation is not seen in the atomic carbon maps nor in the 
CO maps (FiguresEJ. Neither is it seen in the FI R con- 
tinuum at 158 /im observed bv iHowe et alJ l|1991() which 
shows a circular symmetry peaking at IRS 5. This could 
mean that [C 11] does not solely trace PDRs, i.e. the molec- 
ular cloud surfaces. A fraction of [C 11] emission may stem 
from the Hii regions and the diffuse outer halo. In the 
following paragraph, we estimate the contribution by the 
H II regions in W3 Main. 

In addition, the modelled FUV field does not trace 
well the observed distribution of [C 11] . The former shows 
a peak at IRS 5 while the latter shows a north-south ori- 
entated ridge of constant [C 11] emission. This mismatch is 
surprising since [C 11] is thought to b e an accurate tracer 
of the FUV field. iHowe et alJ l)l99l|) showed that a ho- 
mogeneous cloud model cannot reconcile this particular 
discrepancy. Rather a clumpy cloud model has to be in- 
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Fig. 6. Overlay of UV field, [C ii] emission, and 6 cm 
radio continuum. Grey-scale: Iraage of the VLA 6 cm 
radio continuum I Tieftrunk et alJ Il997l) . The box de- 
notes the region observed with KOSMA. Solid contours: 
Contours of the calculated spatial structure of the FUV 
field smoothed to 75" resolution. The positions of the OB 
stars are marked with a +. Dashed contours : extended 
[Cii] emission observed by Howe et al.l Ijl99lh with the 
KAO at 80" resolution. Contour levels range from 40% to 
90% in steps of 10% of the peak intensity of 314 K kms^^. 



voked with dense clumps immersed in a dilute interclump 
medium, allowing the FUV radiation to deeply penetrate 
into the cloud. Here, we try to improve on this analysis 
by using the new observations of [Ci] and CO. 

4.3. [C ll] emission from the ionized medium 

[Cii] emission does partly stem from the Hii regions. 
However, the contribution of [Cii] emission from Hii re- 
gions to the observed [Cii] emission is in most cases 
found to be weak. We have calculated the fraction of 
[Cii] emission from the Hll regions in WSMain using 
the emission measures and electron densities compiled by 
[ Tieftr unk et al. ( 199'3) and the formalism of Russell et al.; 
lll98lh (cf. ICrawford et all 119851: iNikola et al.ll200ll) . We 
found that the contribution is very low, varying between 
0.1 and 0.7%. 

This result appears to be valid also on much 
larger scales. Using the CO BE Milky Way survey, 
IPetuchowski fc BennettI (|l993j) find that H ii regions con- 
tribute only ~ 1 % of the total luminosity of the Galaxy 
in the [C ii] line. 



4.4. ISO/LWS data of IRS5/W3 Main 

We retrieved the FIR ISO/LWS data at IRS 5 from the 
data archive and reanalyzed it as described in Sect . 12. 31 



Table 5. Hii regions and exciting OB stars which were 
used to construct the two-dimensional FUV-field shown 
as drawn contours in FigureEl The sp ectral type (col. 
(2)) is an average from the compilation of lTieftrunk et alJ 
(I99I4). Total luniinositi es (col. (4)) assuming ZAMS stars 
are from |Panagiall)l973|) . Colu mn (5) gives the correspond- 
ing names of infra red sources IjWvnn- Williams et al.ll972i 
iRichardson et al.l ir9891 . 



Hi: 
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Region Type 
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^0 
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(25,21) 


68.0 


IRS 2 


B 


06.5 


(-20,1) 


15.0 


IRS 3a 
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06.5 


(32,-132) 


15.0 




J 


08 


(-87,-129) 


6.5 




D 


08.5 


(-76,34) 


5.4 


IRS 10 


Aa 


08.5 


(18,25) 


5.4 


IRS 2a 


H 


08.5 


(-58,67) 


5.4 




F 


09.5 


(-3,-11) 


3.8 


IRS 7 


C 


09.5 


(-72,27) 


3.8 


IRS 4 


G 


BO 


(-47,0) 


2.5 




E 


BO 


(-37,-18) 


2.5 




Ma-g 


BO 


(0,0) 


7x2.5 


IRS 5 


Ca 


Bl 


(-69,28) 


0.5 


IRS 4 



Sum 



151.0 
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Fig. 7. ISO/LWS spectrum of the IRS 5 position in 
W3 Main The overlapping LWS bands are clearly discern- 
able. The dashed line is the isothermal grey body fit which 
results in Tdust = 53K and /3 = 1.2 for Ay = 40 mag. For 
comparison, the drawn lines show the grey body for 58 K 
and 47 K dust temperature at the same f3 value. The posi- 
tions of several atomic fine structure and high- J CO lines 
are marked. 



4.4.1. Grey body fit 

The far-infrared continuum at IRS 5 is modelled using 
the ISO/LWS data. We assumed an isothermal grey body 
spectrum F\ ~ ri_BA(T'dust)(l — exp~'^'i"=*), where the solid 
angle il was fixed to correspond to the average LWS 
HPBW of - 75", Ba (Tdust) is the Planck function at tem- 
perature Tdust, and the dust opacity is parametrized by 
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Fig. 8. Individual ISO/LWS spectra at IRS 5/W3 Main 
(cf. Fig,[2I), Drawn lines show the fitted instrumental line 
profiles. Integrated line intensities are listed in TableEl 

Table 6. ISO/LWS integrated intensities of the FIR lines 
and the FIR continuum at IRS 5. Integrated intensities 
are listed in units of [10~^ergs~"'^ cm~^ sr~^] and were de- 
rived from the measured flux per beam using the beam 
sizes (col, (4)) and corrections for extended source emis- 
sion given in Table 4,10 in the ISO-LWS Handbook v, 1,2, 
The listed errors only reflect the line fitting uncertainty. 
The intensity of the [Ol](145/xm) line is the average of 
two LWS detectors. The error is in this case the deviation 
from the individual fluxes. 



that all of the FIR luminosity is powered by the 05 star 
associated with IRS 2 (TableEJ , 

Assuming /ftr = 2 x 1,6 10~^(47r)~^G'n erg 
s^^cm^^sr^^ IjKaufman et alJll999() . i,e, assuming that 
the stellar FUV flux is solely responsible for the FIR 
continuum, we estimate the FUV field at IRS 5 to be 
Go = 5,3 10^, This is only 13% of the FUV field estimated 
solely from the spectral types of the OB stars: Gq — 4 10^ 
(Table^ which may again indicate that a part of the stel- 
lar photons escape the cloud without heating the clouds. 
This issue will be discussed again in ChaDter l4,5,4l when 
comparing with the results of the PDR modelling. 

To improve on the analysis of the ISO/LWS contin- 
uum, a better calibration of the LWS bands would be 
needed and the variation of LWS beam sizes (55" — 90") 
would need to be taken into account, 

iQldham et alJ l|l994 » derive slightly higher values of 
2dust = 68 K and f3 = lA from fitting a single-temperature 
grey body spectrum to FIR data at higher angular resolu- 
tions of < 27", indicating again an increase in temperature 
with smaller beam. They also fit an optical extinction of 
Av = 281 mag for a source diameter of 13", Comparing 
this Ay with the Ay we derive within a 75" beam, indi- 
cates a source size of 28", 



4,4,2, Atomic fine structure lines 

Integrated intensities of the transitions detected with 
ISO/LWS at IRS 5 are listed in TableEl 



IRS 5 



A Intcgr,Intensity HPBW 



CO 14-13 
CO 15-14 
CO 16-15 
CO 17-16 



[Cii] 
[Oi] 

[Nil] 
[Oiii] 

[Oi] 

[Niii] 

[Oral 



[^m] 



186, 0,21 ± 0,03 

174, 0,30 ± 0,03 

163, 0,08 ± 0,03 

153, 0,13 ± 0,03 



58, 


1,54 ± 0,03 


45, 


0,77 ± 0,13 


22, 


0,06 ± 0,06 


88, 


3,8 ± 0,2 


63, 


8,3 ± 0,3 


57, 


1,7 ± 0,3 


52, 


19 ± 1 



66, 
66, 
66, 



78, 
76, 
86, 
84, 



FIR-Continuum 



13500 



Tracers of PDRs. The lines of C ii, O i, and CO are promi- 
nent tracers of the PDRs at the surfaces of the molecular 
clouds. Of all FIR lines tracing the neutral medium, the 
63/im line of O i is the major coolant at IRS 5 being a fac- 
tor of 4,1 stronger than the [Cii] line. Moreover, this [Oi] 
transition may be optically thick and thus its intensity 
strongly reduced by foreground absorption. This is indi- 
rectly confirmed by detailed PDR modelling described in 
the Sect, ESI 

The [C ii] integrated intensity dete cted with ISO/LWS 
is 70% of the intensity measured by iHowe et alJ l|l99lh 
with the KAO and thus in agreement within the expected 
calibration error. 



■''dust — Tv{X/Xv) ^. The total II2 column density de- 
rived in Sect 13, 6l corresponds to an optical extinction of 
Ay = 40mag= 1,086 ry ( Bohhn ct al, 1978) and was held 
fixed as well, A least squares fit of the dust temperature 
and the spectral index results in Tdust = 53 ± 1 K and /? — 
1,21 ±0,02 (Fig,[7I), The total brightness, obtained by inte- 
grating over all wavelengths, is 2,1 10^^^ Wcm^^bcam^^ 
or 13,5 erg s~^ cm~^ sr~^. The total luminosity within 
the beam is 3,4 10^ Lq, Comparison with TableEl indicates 



Tracers of the ionized gas. The ions Nil, Oiii, Nlll, 
have ionization potentials larger than 13,6 eV and their 
fine structure lines therefore trace the ionized medium. 
Although, the [Nil] line at 122 /im is not detected, the 
lower limit of the [C 11] /[Nil] line fiu x ratio can be esti- 
mated to at least 27. A ratio of 9 (cf, iNikola et aPbOOlf) 
would be expected for [C 11] being emitted entirely from 
the ionized gas. This indicates that at most 33% of the 
[C 11] emission stems from the diffuse ionized medium. A 
more stringent estimate is derived using the properties of 
the Hii regions in W3Main as discussed in Sect. 14, 31 
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Fig. 9. Integrated CO intensities as a function of upper 
rotational quantum number J for IRS 5 in W3 Main and 
CO cooling curves produced by a non-LTE radiative trans- 
fer model. Rotational transitions upto J ~ 8 were ob- 
served with KOSMA for which we plot the 15% uncer- 
tainty. The high-J CO transitions upto J — 17 were ob- 
served with ISO/LWS. We plot the uncertainties listed in 
TableEl Model 3 is a two-component model adding the 
contributions of the warm gas (model 1) to a cold compo- 
nent at Tkin = 40 K, n = 10'^ cm~^ containing 90% of the 
mass. 

4.4.3. Carbon monoxide emission lines 

Four high-lying rotational lines of gas phase CO were 
detected with ISO/LWS at IRS 5. These lines trace the 
hot and dense gas, upper level energies lie above 500 K 
and critical densities are about lO'' cm~^ and higher. 
To better constrain the physical parameters, we used 
a non-LTE escape probabilit y radiative transfer model 
ijStutzki &: WinnewisseJll985|) assuming a CO abundance 
of 8 1 0~^ and using cross sections from ISchinke et al.l 
l| 19851) . Figureini shows the CO and ^^CO Hue integrated 
intensities observed with KOSMA and ISO/LWS. 

Although the model simplifies the source structure by 
neglecting temperature and density gradients, and vary- 
ing beam filling factors are not considered, common so- 
lutions can be found to fit all 9 observed line intensities. 
The minimum column density to reproduce the observed 
CO fluxes is iV(H2)= 3 lO^i cm^^ only about 10% of the 
total H2 column density at IRS 5 (Table^)). (We assumed 
line widths of 9kms^^ for the ^^CO lines and 6kms^^ 
for the ^■^CO lines.) The column density is largely deter- 
mined by the mid- J CO transitions, rather independent of 



the kinetic temperatures and densities. The kinetic tem- 
peratures need to be high to reproduce the high-J CO 
transitions detected with ISO/LWS while the ISO transi- 
tions are subthermally excited. However, two solutions ex- 
ist which bracket the parameter space of acceptable mod- 
els: 1. a high gas temperature of 200 K and rather low den- 
sity of 3 lO^'cm"'^ and 2. a lower gas temperature of 140 K 
and a higher gas density of 10^ cm~^. Both solutions indi- 
cate a pressure (nT) of (1 - 5) lO'^Kcm'^. The CO 15-14 
transition is only slightly stronger than the CO 7-6 tran- 
sition, the ratio of both intensities is 1.29. On the other 
hand, the intensity ratio of CO 15-14 o ver CO 1-0, de- 
rived fr om the model results, is about SOO. lTieftrunk et al] 
1 1998b^ derived a kinetic temperature of 150 K at IRS 5 



from ammonia observations, consistent with our results. 

However, only about 10% of the gas is at these high 
temperatures and densities. A second phase must ex- 
ist containing the bulk of the gas mass and at signifi- 
cantly lower temperatures. We model this phase assum- 
ing iV(H2)= 31022cm-2, Tkin = 40K, n = lO^cm-^. 
The dashed curve in Figure|Sl shows the combined flux 
(model 3) of the warm phase (model 1) and the cold phase. 
The cold phase hardly contributes to the observed cooling 
fluxes of the mid-J and high-J CO transitions, although 
it contains 90% of the gas mass. 

A [CO]/[^^CO] abundance ratio of 22 is needed to 
fit simultaneously the observed ^'^CO lines. The ^^CO 
8-7 line in particular is crucial in determining the col- 
umn density and abundance ratio. The [CO]/[^^CO] abun- 
dance is significantly lower than the ratio of 89 found by 
iMitchell et alJ Ijl990h at sub-arcsecond resolutions, trac- 
ing the immediate hot and dense environment in front of 
IRS 5. The latter ratio is co nsistent with the canon ical 
interstellar abundance of 65 l|Langer fc Penziaa 1 1990(1 . It 
is unlikely that ^'^CO fractionation plays a major role at 
W3/IRS 5 since the kinetic temperatures are much higher 
than 36 K, the energy gained in substituting ^^C by "'^^C in 
CO. The low [CO]/[^^CO] abundance ratio derived from 
the escape probability model may therefore indicate that 
a more detailed source model is needed, including density 
and kinetic temperature gradients along the line of sight. 
Indeed, in a more realistic model as the PDR model, the 
15-14 and 4-3 CO lines are not produced by gas at the 
same temperature and the same density. 

A similar analysis as we did for IRS 5 combin- 
ing submi llimeter and far-in f rared CO lines was con- 
ducted by iHarris et all ^2^; IStutzki et all (l98?), and 
IStacev et al. (19941) at FIR peaks of the Galactic star 
forming regions M 17, S 106, and in the Orion molecu- 
lar cloud. The derived minimum column densities of the 
warm gas and the kinetic temperatures are very simi- 
lar (7V(C0)= (0.5 - 6)10i*^cm-2, T^n = 100 - 450K, 
(1 - 4)104cm-3) to the solution we find in W3IRS5. 
These conditions appear to be typical for Galactic molecu- 
lar clouds in the vicinity of massive star formation. Similar 
conditions of the warm molecular gas, but on a much 
larger scale of 180 pc are found in the starburst nucleus of 
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NGC 253 by analyzing CO cooling curves l|Bradford et alJ 

l2nn4j) . 



4.4.4. Line cooling from [Cii], [Ci], CO, and [Oi] 



The ratio of [C ii] intensity over the FIR continuum within 
the ISO/LWS beam of ~ 75" at IRS 5 is 1.1 10"'* (cf. 
Tableinj ■ The same ratio averaged over the en tire mapped 
regio n is a factor 14 higher, i.e. 1.5 10~^ (JHowe et al.l 
Il99l[) . Similar values of about 10^^ as found locally at 
W3 IRS 5 have also been found in other Galactic star form- 
in g regions like Orion A , M 17, W 51, or DR 21 (see Table 5 
in lstacevet"ai]ll99l|^ . On the other hand, the CORE 
large scale low resolution map of the Milky Way shows a 
higher and rather const ant ratio with little scatter between 
5010-'' and 17010"'' (JFiiseniraDIliil. [Cii]/FIR ra- 
tio rati os of upto 0.01 have been found in other Galactic 
clouds |Stace}^et al. 1991J) using KAO and in external 
galaxies (jStacev et al...l99ll iMalhotra et al]l200lj) . Based 
on ISO observations. iMalhotra et all l)200lj) found typical 
[Cii] /FIR ratios for normal galaxies (excluding AGNs) in 
the range 0.01 to 0.001. However, there are again signif- 
icant exceptions: two warm and active galaxies show ra- 
tios of less than 2 lO"**, reminicent of the ratio we find at 
IRS 5. A significant [Cii] deficiency has also been found in 
ULIRGs iJLuhman et al.l ll998) . A number of explanations 
are discussed in the literature. One possibility is that in- 
creased Gq/u ratios lead to increased positive grain charge 
which in turn leads to less efficient heati ng of the gas by 
photoelectrons from dust grains. See e.g. iKaufman et al.l 
(.1999.) for a discussion. 
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Fig. 11. Variation of physical parame ters with Ay 
for a PDR model l|Kaufman et alJ Il999l Wolfire priv. 
comm.)with n = 1.810^cm~^ and Go = 10^ which is 
characteristic for IRS 5 (cf. TableO. (a) Local abundances 
n{X)/{n{H) + 2n{H2)) of atomic carbon and CO. The 
CI abundance peaks at Ay ^ 4.2 mag reaching 1.2 10""* 
while the CO abundance rises steadily with Ay reach- 
ing an abundance of 1.410^'' at Ay = 10 mag. (b) The 
gas temperature (drawn line) varies between 4000 K at the 
surface and ~ 70K at Ay — 10 mag. The contours labeled 
CO 15-14 and CI 1-0 show the line intensity integrated 
along Ay from the outside in to greater depths and nor- 
malized to the integrated intensity at Ay — 10 mag. The 
CO 15-14 line stems from regions slightly further out and 
slightly warmer than the region where CI 1-0 is emitted. 
The extinction upto which half of the intensity is emitted, 
is 4.2 mag for CO 15-14 and 4.6 mag for CI 1-0. At these 
extinctions, the gas temperatures are 213 K and 140 K, 
respectively. 



The sum of all observed CO intensities is 1.010 '^erg 
s~^cm~^ sr"-', only a fifth of the total modelled CO line 
cooling for J upto 20 (Model 1). The total CO flux is 
only weakly influenced by self-absorption, since we have 
corrected for this effect to first order, and secondly, the 
CO cooling is mainly produced by the highly-excited lines 
which are less affected by self absorption. The importance 
of the unobserved transitions between J ~ 8 and J = 13 
to the total cooling is evident from FigureEland from the 
escape probability analysis which shows that the peaks 
of the best fitting modelled cooling curves both lie near 
J = 12. 



The total CO cooling by far exceeds the gas cooling 
by the two carbon lines: Aco/^c — 123. The two [Oi] 
lines are more important: Aq/Aco — 1-9, Aq/Acii = 6. 
The relative contributions to the total gas cooling are 
Ago = 31.6%, Aq = 58.7%, Acn = 9.7%, Ac = 0.3%. 
The relative contribution of [O i] may be even stronger 
due to high optical depth and foreground absorption of the 
[O i] (63 /im) line. The total gas cooling fiux is 15.5 lO^'^erg 
s~^cm~^ sr~^ and thus 0.1% of the dust cooling FIR flux. 



4.5. Detailed PDR modelling of the submillimeter and 
far-infrared emission 



We ha ve used the results of PDR models of lKaufman et alJ 
l|1999^ to compare with the observed line ratios. These 
models assume a constant density of a semi-infinite 
plane-parallel slab illuminated by the interstellar ra- 
diation field and improve on the older models of 
iTielens fc Hollenbachl l|l985l) . PAHs are found to domi- 
nate the grain-photoelectric heating and are therefore in- 
cluded. ^'^C is not included in the chemical network. The 
model computes a simultaneous solution for the chemistry, 
radiative transfer, and the thermal balance in PDRs. 

4.5.1. Estimating density and FUV field from the line 
ratios 

Wc have compared the models with selected line fiux ratios 
of the PDR tracers: CO (4-3, 7-6), [Ci] (1-0, 2-1), [Cii] 
158Aim, and, at IRS 5, [OI]{63^nJi, U5^ira) and CO (15- 
14). Figureirni shows the detailed comparison at the four 
selected positions in W3 Main (cf. Fig.EI TablesEl H EJ . 
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Fig. 10. Results from the PDR model of iKaufman et al.1 1 1993) at the four selected positions in W3Main. Contours 
show the observed ratios of line integrated intensities as a function of FUV f lux in units of the Habing field Go and 
the H2 volume density n. We use the [Cii] data from KAO I Howe et al.lll99l|) . The straight horizontal lines show the 
FUV field modelled from the embedded OB stars and smoothed to the resolution of the line data, i.e. 75" (TableEJ. 
Greyscales denote the reduced x^ derived from the observed line ratios, weighted with the uncertainty of 20%, and 
not taking into account the ratios of [Oi](63/im) and CO 15-14 at IRS 5. The cross marks the best solution Xj^, 
Greyscales vary between 1 and 5 times %„ 



/2. . 



The plots span a density range from n = 10 to 10^ cm^'^ 
and a FUV flux range from Gq = 0.3 to 3 10^, i.e. the 
full range of the Kaufman models. Horizontal dashed lines 
mark the beam averaged FUV flux estimated above in 
Sect. l4.l| from the spectral types of the OB stars. 

We do not expect to find one solution of Go,n for 
each observed position, i.e. an intersection of contours at 
one point. This is because, first, the observational cali- 
bration error of 20% of the ratios and, second, the line 
emission surely cannot be entirely described by slabs of 
only one density. The analysis given here will however al- 
low to derive average clump densities, average impinging 
FUV fields, and several other average clump parameters. 
Indeed, the line ratio contours shown in Figure^] for the 
four representative positions do not all intersect. However, 
the observed four ratios of [C ii]/[C i] 1-0, CO 7-6/[C i] 2- 
1, [Ci] 1-0/CO 4-3, and [Cii]/CO 4-3 do almost intersect 
at rici = 2 10^ cm^^ and Go > 10^ at the three positions 
IRS 5, IRS 4, and North. The ratios of CO over [C i] lines 



largely determine the densities and are rather indepen- 
dent of Gq. The high FUV fields are largely determined 
by the [C 11] over [C i] ratios and also by the [C 11] over CO 
4-3 ratios. We find still higher densities of upto 10^ cm^"^ 
from the two ratios of [Ci] 2-1/1-0 and CO 7-6/4-3 at 
all positions. 

At IRS 5, we also obtained the CO 15-14/7-6 ratio. 
Very similar to the CO 7-6/4-3 ratio, it indicates high 
densities of ~ 610^cm~'^ for Go > 10^. However, these 
ratios are still consistent with the best fitting solutions 
within a deviation of better than 50%. FiguretTTl shows 
the variation of abundances, temperatures, and intensities 
with optical extinction for the best fitting PDR model at 
IRS 5: n = l.SlO^cm^^^ Gq = 10^ i.e. Go/n = 5.5 cm^ 
(cf. Tablc|7I). It shows that the bulk of CO 15-14 emission 
stems from a region of ~ 210 K temperature while the 
bulk of CI 1-0 emission stems from a region of ^ 140 K. 
While the detailed PDR modelling shows a strong tem- 
perature gradient, the characteristic CO temperature of 



C. Kramer et al.: Emission of CO, C I, and C ll in W3 Main 



15 



Table 7. Results of detailed PDR modelling at the four selected positions in W3 Main. The minimum y^ was used 
to derive an estimate of the best fitting clump density rid and UV Flux. Tgurf is the surface temperature of the 
corresponding clump predicted by the PDR model. The average density nav, the volume filling factor (^y, the average 
clump diameter Dch and the number of clumps are derived using the PDR model results and the total beam averaged 
H2 column density resulting from the above LTE analysis (TableQ)) . The area filling factor $^ is the ratio of observed 
[C 11] intensities over the modelled intensities. 
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210 K is only slightly higher than the results of the single 
component escape-probability modelling of the mid and 
high- J CO lines described in sec. l4.4.3l (Fig.l^. There, we 
had fitted temperatures and densities of between 140 K at 
10^ cm~^ and 200 K at 310^cm~'^. However, a homoge- 
neous model with 200 K at n = 10*^ cm~'^ would be pre- 
dicting stronger CO lines than considered in Section l¥.4.3l 
Thus, with regard to CO emission, the esc. prob. model 
and the PDR model do not agree with each other com- 
pletely. 

The LTE analysis of the CI line ratio in section lS.4.21 
indicated CI excitation temperatures of > 120 K at IRS 5, 
in rough agreement with the detailled PDR modelling. 

4.5.2. The [Oi](63^m) line is optically thick 

From ISO/LWS at IRS 5, we also obtained ratios of 
[Oi](63/im)/[Cii] and [Oi](145^m)/(63/im). Assuming a 
FUV field of at least Go — 10^, these ratios indicate low 
densities of less than 10"* cm"'^, clearly deviating from the 
above scenario. However, the [O i] ( 63 /im) line is optically 
thick all over the parameter space iJKaufman et al.lll999(l 
and thus probably strongly affected by self-absorption of 
the cold foreground gas component also visible in the low- 
J and mid- J CO lines. The line strengths of the CO 4-3 
and 3-2 lines are reduced by a factor of 1.5 because of self- 
absorption (Fig.|SJ). Decreasing the [Oi](145/im)/(63/im) 
ratio by a factor of 2 would lead to a contour intersect- 
ing with the best fitting Go,n solution of the other ra- 
tios. However, the [Oi](63/xm)/[Cii] would need to be in- 
creased by a factor of 15. Thus, the [Oi](63/^m) line is 
not consistently described in the current model. Velocity 
resolved observations to study the structure of FIR line 
shapes and thus confirm self-ab sorbed profiles will b ecome 
possible with GREAT/SOFIA l|Gusten et al.ll2003(l in the 
near future. 



4.5.3. The quality of the fits 

To quantify the correspondance of the ratios, we have plot- 
ted in FigurelTUl the x^ per degree of freedom, i.e. the re- 
duced x^, of the observed ratios in greyscales. To allow for 



a comparison between the four positions, we did not take 
into account here the ratios of [O i] and CO 15-14 which 
were only obtained at IRS 5. For the six remaining ratios, 
the degree of freedom is 5 — 2. ^ The best fitting solution 
Xmin is marked by a cross and listed in Table[7| A value 
of Xmin '^sar 1 is found at IRS 5 and IRS 4 which signifies 
that the ratios agree well within their observational errors. 
The agreement is worse at the two other positions. 

The best fitting FUV fiuxes of the PDR models arc sur- 
prisingly high, i.e. upto Go = 2 10^, though with a higher 
uncertainty compared to the best fitting densities. At all 
positions, the fitted FUV flux is higher than the beam 
averaged, purely geometrically attenuated FUV fluxes es- 
timated above from the luminosity of the embedded OB 
stars. 

4.5.4. Line and FIR cooling at IRS 5 revisited 

The line fiuxes observed with KOSMA and ISO/LWS at 
IRS 5 had already been discussed in Chapter l4.4.4l Here, 
these fluxes are compared with the absolute fluxes pre- 
dicted by the best fltting PDR model at IRS 5 (cf. TableEI). 
TablelHl shows that model and observations are in agree- 
ment to within 50% for the total fluxes of CO, [Cii], and 
[Ci]. The [Cii] flux observed with the KAO fisted here 
is a factor of 1.5 larger than the flux observed with ISO 
which equals the PDR model result. 

As discussed in Sec. l4.5.!21 the [O i] 63 /im line is pre- 
dicted to be much stronger than observed. The total ob- 
served [O i] cooling is only 10% of the modelled flux. 

As already noted above, the FUV flux derived by fit- 
ting the line ratios to the PDR model is Go — 10^, surpris- 
ingly high relative to the beam averaged FUV fiux derived 
from the OB stars. Go = 4 10^, and higher than the FUV 
fiux derived from the beam averaged FIR field. Go = 5 10"* 
(Sec. 14. 4.1)) . Correspondingly, the FIR fiux derived from 
the modelled Go via Jftr = 2 x 1.6 10-3(47r)-iGo erg 
s~^cm~^sr~^ (Kaufm an et alJ[l999,) . is high, /fir = 260 



^ Only 5 ratios are independent: the ratio CO 7-6/CI 2-1 can 
be derived from CI 2-1/CI 1-0, CI 1-0/CO 4-3, CO 4-3/CO 7-6. 
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Table 8. Observed fluxes at IRS 5 in units of 10 '^ erg 
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sr^ compared to the best fitting PDR model 
', Go = 10^. The total gas cooling is 
gas — ^icii -t- ^i-ci -r Aco + A-oi- " The total "observed" 
CO cooling Aco ^^-s derived from the esc. prob. model 
fitted to the observed fluxes (cf. Chapter l4.4.4|l . 
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ergs ^ cm ^ sr ^, i.e. a factor 20 higher than the observed 
flux. 

This indicates that the line emission is mainly excited 
in the dense inner regions near the embedded OB stars 
where the FUV flux is higher. The FIR continuum on the 
other hand also stems from the colder outer surroundings 
thus leading to a lower average FIR intensity. 



4.5.5. A clumpy PDR model 

Using these results in combination with the total beam av- 
eraged H2 column density resulting from the LTE analysis 
at the four positions (TableEJ , allows to derive more phys- 
ical parameters characterizing the line emitting gas. We 
assume here, that line emission arises in individual clumps 
immersed in an interclump medium which is not emitting 
in the observed transitions. Motivated by the morphol- 
ogy of the KOSMA maps shown in Figure|21 we assume 
conservatively that the line of sight extent of the cloud is 
.Dcioud = 0.84 pc (75"). This allows to use the total column 
density to derive a rough estimate of the average volume 
density Uav per beam. The ratio of Uav with the average 
clump density rid derived above, gives the volume filling 
factor of the clumps <I>y. The area filling factor $^ of the 
clumps can be estimated from the ratio of observed inten- 
sities Jobs over the modelled intensities /mod - it can be 
larger than unity for optically thin emission. The clump 
diameter Dd can be estimated from Ddoud^v/^A- The 
number of clumps is then given by <I'v'(-Dcioud/£'ci)^- 

Tabled summarizes the resulting quantities at the four 
positions. Results are much less certain for the southern 
position than for IRS 5. Only less than 10% of the vol- 
ume is filled with a few hundred clumps per beam area. 
Clumps have a typical diameter of 0.05 pc, i.e. 4.4", and 
are thus unresolved. The resulting typical clump mass is 
M = rici 4/3 vr {Da/2f = 0.44 Mq. 

The area filling factor derived from [Cii] is ~ 1 ± 0.5. 
Area filling factors and clump diameters are not signif- 
icantly changed when e.g. using the modelled [Ci] 1-0 
intensities instead. Another reason for discrepancies be- 
tween the mode lled and the obs e rved intensities are geo- 
metrical effects. Kaufman et al.l l|1999|) assume a face-on 
geometry, a viewing angle of < 90° would lead to an in- 
crease in the intensity of optically thin lines. 



We cannot exclude the possibility that a fraction of the 
CI and CII emission stems from the interclump medium 
since the CI and CII transit ions have critical deii sities of 
only 5 10^ to 3 10^ cm'^ (cf. iKaufman et"al]ll999l) 

These results confirm the results of lHowe et all lll99lh 
which had been based solely on the [C 11] map and on 
the estimated FUV field. They had found that a cloud 
clumpy model can be adjusted to be roughly consis- 
tent with the observed [Cil] morphology. Their best fit- 
ting model is very similar to ours: their clumps have a 
density of rid ^ 10^ cm""^, the volume filling factor is 
$y -0.05-0.07. 

Further indirect support for our clumpy source model 
comes from the analysis of CS 2-1, 5-4, and 7 6 FCRAO 
and K OSMA small maps of WSMain bv lOssenkopf et all 
1)200 1|) . The observed line profiles are best fitted by an av- 
erage density of 0.7 10'* cm~'^ within a radius of 1 pc while 
the minimum central density reaches 510^cm~'^, higher 
than the clump density we have derived from the PDR 
analysis, but similar to the high densities inferred from 
the [Ci] 2-1/1-0 and CO 7-6/4-3 ratios. 

Direct evidence for the existence of very dense and 
small molecular clumps is delivered by VLA maps of 
Ammonia in W3 Main at 3" resolution bv lTieftrunk et alJ 
l|1998af) . These show several clumps of ~ 0.02 pc radius 
and densities of a few 10^ cm~^ in the vicinity of IRS 5 
and IRS 4. 

Still higher densities at temperatures of almost 600 K 
are known to exist on scales of less than 0.01 pc (< 
2000 AU, < 1") in the immediate vicinity o f IRS 5 from 
the absorption study of iMitchell et al.1 (|l99ff). 

Our scenario of many small, high density clumps with 
a small volume, but a high area filling factor should never- 
theless not be taken too literal. It simply implies that the 
cloud structure is very far off a homogeneous cloud, but 
shows a very much broken up structure with high density 
contrasts and steep density gradients within the cloud. 



5. Summary 

1. We have mapped the W3Main region in the two fine 
structure lines of atomic carbon and in several mid- J 
CO transitions using KOSMA. All maps show a peak 
of integrated intensities near IRS 5 and IRS 4 with a 
steady drop of intensities in all directions. Spectral 
line maps show structure down to the spatial resolu- 
tion limit of > 42", i.e. > 0.5 pc. To obtain a common 
spatial resolution, all data discussed below are at a 
common resolution of 75" (0.84 pc). In addition, we 
have obtained ISO/LWS data at IRS 5. 

2. We find a constant CO 7-6/4-3 line integrated inten- 
sity ratio of — 1 at all positions in W3 Main after cor- 
recting for self-absorption seen in the *^C0 lines in the 
vicinity IRS 5. Assuming thermalized emission, this ra- 
tio indicates excitation temperatures in excess of 85 K 
and densities in excess of — 10^ cm~^ for the gas emit- 
ting these mid-J CO lines. A more detailed non-LTE 
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analysis however reveals that densities of ~ 10^ cm~'^ 
are sufficient to explain the observed CO cooling fluxes. 

3. The observed high [Ci] 2-1/1-0 line ratios of between 
1.3 and 2.9 indicate optically thin carbon emission and 
excitation temperatures of about 100 K. 

4. LTE analysis of carbon column densities show that the 
C i/CO abundance ratios stays constant at ~ 0.11. The 
fractional abundances of C Ii/C i/CO are ~ 20 : 10 : 
70 at IRS 5, rising to ^ 32 : 7 : 61 in the south of 
W3 Main. 

5. A model of the FUV field distribution, derived from 
the embedded OB stars, assuming pure geometrical 
attenuation, ignoring the cloud density structure, and 
smoothed to 75" resolution, shows a peak field at IRS 5 
of Go = 10^'^ dropping to values below 10^'^ in the 
outskirts of W3 Main. 

6. The north-south orienation of the [Cii] emission 
l|Howe et al.lll99ll) is not seen in the maps of [C i] and 
CO emission. Nor is it seen in the modelled FUV field 
distribution. 

7. At IRS 5, the (average) dust temperature of the FIR 
emission detected by ISO/LWS is 53 K with a spec- 
tral index oi (3 — 1.2. The [Oi] 63 /im line is a factor 
4 stronger than the [Cii] 158 ^m line. The [Cii]/FIR 
flux ratio is 0.01%. We estimate that less than 1% 
of the [Cii] emission stems from the Hii regions in 
W3Main. The relative contributions to the total gas 
cooling are Aco = 31.6%, Aq = 58.7%, Acii = 9.7%, 
Ac = 0.3%. Aco refers to the total CO cooling mod- 
elled with the escape probability radiative transfer 
model in Sec. l4.4.4l All other numbers refer to the ob- 



served values. The ratio of A, 



ras/AHiist, 1 

of lKaufmg 



is 0.1%. 



8. We have used the PDR model of lKaufman et aO l)l999|) 
to estimate the average densities of the emitting region 
and the impinging FUV fields from the observed line 
ratios at four selected positions. The [O i](63/im) line is 
probably optically thick and strongly affected by self- 
absorption in the colder foreground material at IRS 5, 
and thus discarded from the further analysis. The re- 
maining six line ratios indicate a common solution of 
Go ^ 10^ and rici ^ 2 10^ cm"'^ at all four positions. 
The reduced x^ is low at IRS 5, IRS 4, and the north- 
ern position, but much higher at the Southern position 
where no well defined common solution is found. We 
have assumed that all line emission arises in individual 
clumps immersed in an interclump medium which is 
not emitting in the observed transitions. We find that 
the emitting gas consists of a few hundred clumps per 
0.84 pc beam, filling between 3 and 9% of the volume, 
with a typical clump radius of 0.025 pc (2.2"), and typ- 
ical mass of 0.44 Mq. 

9. The line ratios of [Ci] 2-1/1-0 and CO 7-6/4-3 are 
consistent with very high clump densities of upto 
10^ cm""^ in the framework of the PDR model. This 
finding may be explained by a more detailed model of 
clumps having a variation of densities (and masses). A 
more detailed model of W3 Main would also need to 
take into account absorption by foreground gas. 



Although the discussed model scenario is not able to 
reproduce each detail of the line intensities and their vari- 
ation within the cloud, it is remarkable that the scenario 
of a UV penetrated clumpy cloud is able to reproduce the 
relative (and to a large degree also the absolute) intensities 
of some dozen lines simultaneously with only a relatively 
small set of model parameters such as the clump average 
density and size, most other parameters being constrained 
by the UV intensity distribution and the temperature and 
chemical gradients calculated through the PDR models. 
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